ITEP-LAT/2007-08 



Center vortices as rigid strings 

P. V. Buividovich* 
Belarusian State University, 220080 Belarus, Minsk, Nezalezhnasti av. 4 

M. I. Polikarpoyt 
ITEP, 117218 Russia, Moscow, B. Cheremushkinskaya str. 25 
(Dated: May 25, 2007) 

It is shown that the action associated with center vortices in SU (2) lattice gauge theory is strongly 
correlated with extrinsic and internal curvatures of the vortex surface and that this correlation 
persists in the continuum limit. Thus a good approximation for the effective vortex action is the 
action of rigid strings, which can reproduce some of the observed geometric properties of center 
vortices. It is conjectured that rigidity may be induced by some fields localized on vortices, and a 
model-independent test of localization is performed. Monopoles detected in the Abelian projection 
are discussed as natural candidates for such two-dimensional fields. 
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I. INTRODUCTION 

Yang-Mills theory is often believed to be equivalent to 
some string theory, however, up to now there is no way to 
detect thin strings behind physical chromoclectric string 
of finite thickness, which gives rise to linear potential be- 
tween test colour charges and which is usually seen in 
numerical simulations [1, 2]. On the other hand, closed 
magnetic strings, or center vortices, can be directly de- 
tected [3] and seem to be thin [4]. A model-independent 
argument in favor of physically thin center vortices in 
continuum pure Yang-Mills theory was proposed recently 
in [.'i]. It is known that the full QCD string tension is 
reproduced with sufficiently good precision if one consid- 
ers only the contribution due to topologically nontrivial 
winding of center vortices and Wilson loop [-i, (i]. More- 
over, in [3, 4] it was demonstrated that the area of center 
vortices in SU (2) LGT scales in physical units of length. 
These facts imply that center vortices are not just lat- 
tice artifacts, but rather correspond to some physically 
significant objects. Center vortices arc usually detected 
using the maximal center gauge and arc seen as closed 
self-avoiding surfaces, which occupy only a small fraction 
of lattice plaquettes [3, 4]. Typically one observes a large 
percolating vortex, which extends through the whole lat- 
tice, and a number of small satellite vortices, whose size 
typically does not exceed several lattice spacings [7] . 

As the total physical area and size of center vortices re- 
main finite in the continuum limit, there should be also a 
continuous description of vortex geometry, which is char- 
acterized by some finite Hausdorf dimension ■ Perco- 
lating vortex by definition has Hausdorf dimension equal 
to 4. In order to define the dimensionality of small satel- 
lite vortices in SU (2) LGT, their size L was measured as 
the function of their area S. The size of the vortex was 
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defined as the maximal distance between points which 
belong to the vortex. Average size of small vortices as 
the function of their area (in lattice units) is plotted on 
Fig. 1. Fit of the form L = const • S^^'^ (solid curve 
on Fig. 1) gives d = 1.9 ± 0.1. For small vortex areas 
{S < 30) this number simply refiects self-avoiding of vor- 
tex surfaces, but for larger values of area this fit indicates 
that small vortices tend to be smooth surfaces. 
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FIG. 1: Average size of center vortices as the function of their 
area {f3 = 2.60, 28^* lattice) 

How can one describe the properties of such surfaces? 
A necessary condition for physical scaling of area of ran- 
dom surfaces is the cancelation between their entropy, 
which for self-avoiding surfaces grows linearly with the 
area of surface in lattice units, and the bare string ten- 
sion, which therefore should be finite in lattice units [(S]. 
A remarkable result of [4] is that the excess of action 
associated with center vortices is indeed proportional to 
the of vortices in lattice units. However, for the 
simplest model of random surfaces with Nambu-Goto ac- 
tion such naive balance between action and entropy does 
not lead to physical surfaces because of branched poly- 
mer problem [s] , therefore one has to consider some more 
complicated model. 

A well-studied model which can probably describe 
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smooth surfaces even in four and three dimensions is the 
model of rigid strings [8, 9, 10, 11, 12]. The model was 
first analyzed in [9, 10], where it was shown that depend- 
ing on the form of the /3-function of the model either 
the branched polymer phase or the phase of smooth sur- 
faces can be observed. Namely, if the /3-function has no 
nontrivial fixed points, the model reduces to the usual 
Nambu-Goto string flawed by branched polymer prob- 
lem, but if nontrivial fixed point exists, in the vicin- 
ity of this point the model should describe smooth sur- 
faces with finite Hausdorf dimension [9] . Numerical sim- 
ulations indeed confirmed the existence of the phase of 
smooth surfaces in this model [8, 12, 13]. Nontrivial scal- 
ing in the vicinity of the phase transition was observed in 
[12], but in more recent simulations of a related model of 
tethered surfaces first-order phase transition was found 
[13]. However, even if the model of rigid strings has no 
continuum limit, it can still serve as an effective model, 
for which UV cutoff is set by some other fields in the the- 
ory. It is known, for instance, that rigidity terms arise 
in the effective string action after integrating over world- 
sheet fermions [S, 14, 15] or some four-dimensional mas- 
sive fields coupled to the string worldsheet [11, IG, 17]. 

The aim of this paper is to fit observed properties of 
center vortices in the model of rigid strings. This model 
was first conjectured to describe center vortices in [4, 18]. 
In this paper the correlation between the excess of action 
associated with vortices and the geometry of vortex sur- 
faces will be studied systematically. In the section II it 
will be shown that the effective vortex action can be in- 
deed approximated by the action of rigid strings and that 
the coupling constants before extrinsic and internal cur- 
vatures are finite in the continuum limit. In the section 
III this action is discussed as an effective action induced 
by some two-dimensional fields localized on the surfaces 
of vortices and a model-independent check of localiza- 
tion is performed. Monopoles in the Abelian projection 
of the theory [4, 19, 20] are then discussed as appropri- 
ate candidates for such localized fields in the section IV. 
The dependence of the bare string tension on the lattice 
spacing can be explained if one takes into account the 
contribution of percolating monopole cluster. Interac- 
tion between monopoles, which can be approximated by 
the Yukawa potential with physical mass [21], can also 
partially account for surface rigidity [11, 16, 17]. Finally, 
implications and possible extensions of obtained results 
are discussed. 



11. EFFECTIVE ACTION AS THE FUNCTION 
OF VORTEX GEOMETRY 

In order to measure the correlation between geometric 
properties of vortex surface and the action density cen- 
ter vortices were detected in SU (2) LGT by imposing the 
direct maximal central gauge (DMC) [4]. Simulated an- 
nealing procedure was used to locate true minima. In this 
work the same set of lattice configurations as in [4] was 




FIG. 2: Average excess of action per site as the function of 
extrinsic curvature and lattice spacing 




FIG. 3: Average excess of action per site as the function of 
internal curvature and lattice spacing 



used. The lattices used were of the size 28* for /3 = 2.60 
and (3 = 2.55, 24* for f3 = 2.50, /3 = 2.45 and /? = 2.40 
and for /3 = 2.35. Lattice spacing was fixed by setting 
the value of QCD string tension to ^/a = 440 MeV. 

Local geometry of vortex surfaces was characterized by 
the two simplest local invariants - internal and extrinsic 
curvatures. Internal curvature in lattice units for hyper- 
cubic lattice was defined as a^Rs = 4,—ns, where is the 
number of neighbors of the site s and a is the lattice spac- 
ing [n]. Extrinsic curvature for smooth surfaces can be 
written as K = Ax^Ax^^, where A/ = -^da {\/gg'^''dbf) 
is the two-dimensional Laplacian on the surface and 
gab = daX^dbx'^ IS the induced metric on the surface. 
In order to define extrinsic curvature on hypercubic lat- 
tice two-dimensional lattice Laplacian on the surface was 
defined as a^Afs ~ ngfs — X^/s'j where s' are the sites 

s' 

adjacent to s. Extrinsic curvature in lattice units is then 
a'^Kg = Ax^Ax'^ [8]. The excess of action associated 
with some lattice site on the vortex was defined by aver- 
aging the excess of action over all vortex plaquettes which 
contain this site. The excess of action on plaquette was 
defined as ASp = /3 (l - ^Tr Up) - (/? (l - ^Tr Up)) = 
f ((Tr Up) — Tr Up). Average excess of action per site as 
the function of lattice spacing and extrinsic and internal 
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FIG. 4: Excess of action per site as the function of extrinsic 
curvature at different lattice spacings. Solid lines are plotted 
using the value at a^Ks = 3 averaged between a^Ks = 2 and 



curvatures (in lattice units) is plotted on Figs. 2 and 3 
respectively. 

It can be seen that the excess of action increases with 
both extrinsic and internal curvature, therefore the ac- 
tion should in general depend on both internal and ex- 
trinsic curvatures. Standard dimensional analysis shows 
that only the terms linear in extrinsic and internal cur- 
vature are relevant in the continuum limit. It should be 
also noted that because center vortices are defined on 
hypercubic lattice, typical values of curvature diverge as 
in the continuum limit {a?Ks and a^Rg are finite 
integer numbers, therefore Kg ^ and Rg ^ a^^). A 
simple estimation shows that the number of points on 
the surface with K ^ or R scales to zero as a^, 
which ensures the smoothness of surfaces and the finite- 
ness of the total contribution of bent plaquettes to the 
physical action. Finite values of Kg or Rg for physical 
surfaces near the continuum limit can be in principle de- 
fined by averaging the curvature over physically small 
regions, whose area is nevertheless very large in lattice 
units. 

A peculiar feature of the dependence of action excess 
on extrinsic curvature is the strong peak near a^Kg = 3. 
Such value of extrinsic curvature corresponds to an edge 
where three plaquettes join. It was found that when this 
peak is neglected, the resulting function appears to be 
almost linear in a^Kg (Fig. 4). Presumably the peak 
at a^Kg = 3 is a lattice artifact and only reflects the 
anisotropy of hypercubic lattice. In order to obtain bet- 
ter fits the excess of action at a^Kg =3 was replaced by 
the average between a^Kg = 2 and a^Kg =4 (solid line 
on Fig. 4), which yields almost linear dependence. Bare 
string tension in lattice units and the coefRcient before 
extrinsic curvature were obtained as the intercept and the 
slope of the linear fits of the data on Fig. 4. One could in 
principle try to fit the excess of action per plaquette by 
some polynomial in a^Kg, but anyway dimensional analy- 
sis shows that only the term linear in a^Kg would survive 
in the continuum limit. Nevertheless, the fitting method 



affects numerical values of the coefRcient k (a) before ex- 
trinsic curvature. In general, increasing the number of 
degrees of freedom in fitting functions changes the value 
and the uncertainty of this coefficient, however, these 
values agree within error range when extrapolated to the 
continuum limit. For instance, in order to check the sta- 
bility of fits the data plotted on Fig. 4 were fitted by 
first and second-order polynomials in a? Kg. The values 
of bare string tension (Tq («) (constant term in the fits) 
obtained from both fits agree with a very good precision. 
At finite lattice spacings discrepancies in the values of 
K, are larger, but extrapolation to the continuum limit 
gives consistent values n = 0.065 ± 0.006 (linear fit) and 
K = 0.048 ± 0.014 (quadratic fit). The coefficient before 

{o?Kg^^ contains very large errors and is close to zero. 
This coefficient can only be important in the continuum 
limit if it contains divergences of order a~^, which is not 
likely. 

In order to extract the term linear in internal curvature 
from the data plotted on Fig. 3 the excess of action was 
fitted by a third-order polynomial in a^Rg = 4 — rig. The 
coefficient before a?Rs in this polynomial was assumed to 
be the coefficient before internal curvature in the physical 
action. 

Finally, after omitting all terms which become irrele- 
vant in the continuum limit, for sufficiently small lattice 
spacings a one can write the action associated with center 
vortices in the following form: 

W[S]^ J d^^^{ao{a)Aljy+j{a)R + K{a)K) (1) 
s 

where ^/g = V det gat is the invariant element of area 
on the surface, gab = is the induced metric on 

the surface and Ajjv ~ a^^ is the UV cutoff scale of the 
theory. 

The coefficients ctq (a)i 7 (a) and k (a) as the functions 
of lattice spacing are plotted on Figs. 5, 7 and 6 re- 
spectively. Extrapolation to continuum limit gives the 
following values: 

cro(O) =0.192 ±0.006 
k(0) = 0.066 ±0.003 

7(0) = 0.08 ±0.02 (2) 

Thus coefficients k and 7 are finite in the continuum 
limit and therefore the dependence on internal and ex- 
trinsic curvatures is physical. Quadratic divergence in 
the bare string tension is crucial for the existence of 
smooth physical surfaces, as explained above, and should 
be compensated by a similar divergence in the entropy 
of random surfaces. It is interesting to note that the 
value of bare string tension in lattice units (2), which is 
obtained after taking curvature into account, is smaller 
than the value of action excess (A5p) = 0.540 ± 0.004 
obtained in [4]. The fact that after proper redistribu- 
tion of action excess among operators with appropriate 
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dimensions the string tension is strongly decreased in- 
dicates that the terms with extrinsic curvature play an 
important role in the dynamics of center vortices. 

The action (1) corresponds to the model of rigid strings 
[8, 9, 10]. While the surface entropy is canceled by the 
divergent bare string tension, branched polymer problem 
is circumvented due to the third term. It is interesting to 
note that as the consequence of the Gauss-Bonnet iden- 
tity / (f^y/gR = 2ttx [S] (x [S] is the Euler characteristic 
s 

of the surface), the coefficient before internal curvature 
is proportional to the string coupling constant, which is 
therefore also finite in the continuum limit. 
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FIG. 5: Bare string tension in lattice units tJo (a) 
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FIG. 7: Coupling constant before internal curvature 7 (a) 



directly observed in lattice simulations - for instance, it 
was found that in the maximal Abelian gauge almost all 
trajectories of Abelian monopoles belong to center vor- 
tices [20] . There also exist classical string solutions where 
monopoles are localized on string worldsheets [24, 2.^i]. In 
[26, 27] it was shown numerically that eigenfunctions of 
covariant four-dimensional Laplace and Dirac operators 
are also localized near center vortices. 

If only two-dimensional fields and their interactions are 
responsible for the excess of action on the vortices, total 
excess of action on the vortex W [S] should be treated as 
the effective string action Weff [S] obtained by integrat- 
ing over the fields at fixed geometry: 




exp(-I^e// [S])^ / P0exp 
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FIG. 6: Coupling constant before extrinsic curvature k (a) 



III. FIELDS LOCALIZED ON CENTER 
VORTICES 

From the point of view of field theory quadratically 
divergent term J (P^^(Tq (a) Afjy in the vortex effec- 

s 

five action indicates that some fields are localized on 
center vortices and become effectively two-dimensional. 
Quadratic divergence in the effective action (1) corre- 
sponds to vacuum oscillations of these two-dimensional 
fields [22, 23]. Localization on center vortices was also 



d^^VgLi^] (3) 



where Vcj) is the covariantly defined path integral mea- 
sure and L [(/>] is the Lagrangian of the field (j). Such 
effective action will necessarily contain quadratic diver- 
gence of the form / d'^^y/g Afjy due to vacuum oscilla- 
tions of these two-dimensional fields plus some geometry- 
dependent terms. For instance, if the fields (j) are four- 
dimensional Dirac fermions, the effective action (3) de- 
pends on the extrinsic curvature of the surface [8, 14, 15]. 

In order to check whether some two-dimensional fields 
propagate along the vortex or not, one can consider the 
points on the vortex which are very distant in terms of 
internal geometry on the vortex but are very close in four- 
dimensional space. If the relevant fields propagate only 
along the vortex, correlation between excess of action in 
such points should be much less than between plaquettes 
separated by the same distance along the vortex. In order 
to check this numerically center vortices were represented 
as graphs, each node of the graph corresponding to some 
lattice plaquette. The correlation between action den- 
sities was measured for the points which are separated 
by only one lattice spacing in four-dimensional space 
but no less than 6 spacings along the vortex {d4 < 2, 
(i2 > 6). The standard breadth first search algorithm for 
unweighted graphs was used to measure the distances on 
vortex. In order to reduce anisotropy the sites which be- 
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long to the vortex were linked not only along lattice links, 
but also along the diagonals of plaquettes. As the dis- 
tance measured was used only for lower-bound estimates 
of distances, there was no necessity to use much more 
time-consuming search algorithms for weighted graphs 
such as Dijkstra algorithm. For comparison the corre- 
lation between neighboring vortex plaquettes (^4 < 2, 
d2 < 2) was also measured. Correlation between plaque- 
tte variables Tr Up and Tr Up' was characterized by the 
correlation coefficient p [Tr Up, Tr Up']: 



p [Tr Up, Tr Up' 



{Ti Up Tr Up') -{Ty Up 



((Tr Up) 



(Tr Up)^ 



(4) 
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FIG. 8: Correlation between neighbouring plaquettes with 
d4 < 2, d2 > 6 and d4 < 2, d2 < 2 

The results of these measurements are plotted on Fig. 
8. It can be seen that the correlation in four-dimensional 
space is notably smaller than along the surface of the vor- 
tex, and therefore the fields which are responsible for the 
excess of action on vortices arc more likely to propagate 
along their surfaces. Unfortunately due to insufficient 
statistics it was not possible to measure the correlation 
lengths which correspond to propagation along the vor- 
tex and in four-dimensional space. The latter can be es- 
timated as the inverse mass of the lowest glueball (^ 1.5 
GeV), which is comparable with the inverse lattice spac- 
ing. A measurement which is somewhat similar to the 
one described above was performed in [4], where the av- 
erage excess of action on plaquettes very close to center 
vortices was shown to be zero. 



IV. ABELIAN MONOPOLES AND CENTER 
VORTICES 

Up to now the only way to see directly the content of 
the conjectured two-dimensional field theory is to impose 
maximal Abelian gauge and to look at the trajectories 
of Abelian monopoles [19, 22, 23, 28], which populate 
densely the surfaces of center vortices. In real simula- 
tions about 95% of monopole trajectories belong to cen- 
ter vortices [19, 20]. It is natural to conjecture that two- 



dimensional field theory living on vortex surfaces should 
describe these monopoles upon first quantization. 

Another nontrivial fact which supports the statement 
about the role of monopoles in the dynamics of center 
vortices is the dependence of bare string tension in lattice 
units on lattice spacing (see Fig. 5), which can be approx- 
imated by a linear function ctq (a) = A + Ba with good 
precision. Besides quadratic divergence, this also yields 
a^^ ~ Auv divergence in the bare string tension in phys- 



ical units: 



CTo (a) 



+ a ^B, where A is given 



by (2) and B — 2.2 ± 0.1/m^^. Such divergence usu- 
ally corresponds to self-energy of one-dimensional objects 
[4, 8, 22, 23]. If the density of one-dimensional objects 
per unit of vortex area piu scales in physical units and is 
constant in the continuum limit, while the bare mass of 
these objects is UV divergent and is close to the critical 
value rribare ~ cT^ In 4 [8], for p\D one then obtains the 
following estimation pxu = = 1.5(8)/m~^. Taking 
into account that the density of vortices is py = 24/to^^, 
it is easy to estimate the density of one-dimensional ob- 
jects in four-dimensional space as p'ljj = 37.(9)/m~'^, 
which is in a good agreement with the density of perco- 
lating Abelian monopoles pm = 31.(l)/m~'^ [2s]. The 
difference may be explained by incomplete detection of 
Abelian monopoles and by curved geometry of vortex 
surfaces. 

Geometric properties of monopole trajectories in 
SU{2) LGT were studied in [28]. It was found that 
the properties of monopole trajectories at hadronic 
scale are not described by random walks with dn = 2, 
as in the case of free scalar particles, but rather by 
smooth random walks. Namely, the correlation between 
the directions of tangent vectors to monopole paths is 
characterized by a correlation length w 300 McV, 
which remains finite in the continuum limit. Smooth 
random walks are in the same universality class as the 
trajectories of spinning particles [8], therefore the fields 
associated with monopoles living on center vortices 
are more likely to have nonzero spin. The simplest 
physical model which leads to smooth random walks is 
the propagation of massive Dirac fcrmions in Euclidean 
space [8, 29]. If the monopoles are assumed to be 
Dirac fermions, the mass of the fermions can be roughly 
estimated as ~ 1.5 GeV from the measurements of 
monopole current correlators in [28]. A general smooth 
random walk corresponds to random walk of the tangent 
vector on the three-sphere and therefore includes 
components of all spins, but it is not clear how such 
model can be incorporated into the models of random 
surfaces. On the other hand, random surfaces with four- 
dimcnsional fermions are well studied. It is known that 
the effective action of fcrmionic strings includes rigidity 
terms [14, 15, 19]. As in the model of fcrmionic strings 
worldsheet fermions are massless Dirac fermions, it is 
natural to conjecture that the effective action induced by 
massive worldsheet fermions also induces string rigidity. 
In general, one can expect that for massive fermions 
the effective action besides the local terms of the form 
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(1) should also contain nonlocal terms of the form 

/ d^eivVte) / rf'6\/5(6)i^ (Ci) (fi, 6) K (6), 
where A ^ (^1,^2) is the kernel of the inverse Laplace 
operator on the surface of the vortex. Unfortunately, 
at presently available lattices it is very difficult to trace 
such terms in the effective action. 

The effective Lagrangian governing the dynamics of 
monopoles was obtained in [21] using the inverse Monte- 
Carlo method. It was found that the effective monopole 
action, besides the usual kinetic term, contains Yukawa 
interaction with physical mass as well as four-point and 
six-point interactions. Higher-order interaction terms 
were found to be very small. Taking these results to- 
gether, a reasonable conjecture is that at hadronic scale 
monopoles behave as massive Dirac fermions living on 
center vortices. Yukawa-type interaction may be induced 
if these fermions are coupled to some four-dimensional 
massive field. Coupling to massive four-dimensional 
fields also leads to rigidity terms in the effective action 
[11, 16, 17], and thus such interaction can be also ab- 
sorbed in the effective action (1). 

V. CONCLUSIONS 

In this paper the relation between local vortex geome- 
try and the action density was studied. Direct measure- 
ments of action imply that a reasonable approximation 
for the effective action of center vortices is the action of 
rigid strings [9, lo]. This action can reproduce observed 
smoothness of vortex surfaces and presumably the physi- 
cal scaling of vortex area. The latter possibility depends 
crucially on the form of nonperturbative /^-function of 
the model of rigid strings and at present time can only 
be checked numerically [12]. Unfortunately the values of 
the parameters <jq, k and 7 (2), obtained by extrapola- 
tion to the continuum limit can not be compared with 
the corresponding critical values obtained from indepen- 
dent simulations, because most numerical investigations 
of the model of rigid strings dealt with three-dimensional 
case [12, 13]. As the existence of the continuum limit of 
the model of rigid strings has not been proven exactly, 
it is not clear whether the effective vortex action (1) can 
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